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We show that finite Fermi systems governed by a mean field and a few-body interaction generically 
possess spectral fluctuations of the Wigner-Dyson type and are, thus, chaotic. Our argument is 
based on an analogy to the metal-insulator transition. We construct a sparse random-matrix 
ensemble ScE that mimics that transition. Our claim then follows from the fact that the generic 
random-matrix ensemble modeling a fermionic interacting many-body system is much less sparse 
than ScE. 



I. INTRODUCTION 

Finite fermionic many-body systems (atoms, 
molecules, nuclei) often display spectral fluctuation 
properties that agree with predictions of random-matrix 
theory (RMT), more precisely, with those of the Gaus- 
sian Orthogonal Ensemble (GOE) This fact is 
commonly taken as evidence for chaotic motion 0, 
In RMT, all pairs of states are coupled by independent 
matrix elements. In a many-body system that situation 
arises only in the presence of many-body interactions. 
But atoms, molecules and nuclei are governed by a mean 
field with residual interactions that are predominantly 
of two-body nature. Therefore, an important question 
is: Does a two-body interaction (or, more generally, 
a fc-body interaction with k > 2 integer but smaller 
than the number m of fermions) generically give rise 
to chaotic motion? The question has received much 
attention (see the review j4j]). In condensed-matter 
physics it has recently been addressed as the problem of 
many-body localization p|. 

While for small values of m the question can be an- 
swered by matrix diagonalization, a general answer is dif- 
ficult to obtain both analytically and numerically because 
the Hamiltonian matrices become ever more sparse with 
increasing to: In every row and column the ratio of the 
number of non-zero matrix elements to the total num- 
ber of matrix elements tends asymptotically (to — > oo) 
to zero. Therefore, analytical arguments like the ones in 
Ref. @ (where that ratio is taken to be asymptotically 
finite) do not apply. And for the physically interesting 
cases k > 2 matrices with dimensions that are numeri- 
cally practicable are very far from the sparse limit. 

The most determined effort so far to overcome these 
difficulties was made in Ref. 0. The authors consid- 
ered a fermionic many-body system governed by random 
one- and two-body interactions. For fixed parameter 
values of the one-body part, they determined the critical 
strength of the two-body interaction where the crossover 
from Poisson statistics to Wigner-Dyson statistics takes 
place. Their arguments are based on a combination of a 
perturbative approach and numerical results. Their re- 



sult implies that a system governed by a pure two-body 
interaction is chaotic. Needless to say, the matrices stud- 
ied numerically were far from the sparse limit. 

In this paper we aim at a further clarification of the 
issue. We consider the matrix representation of a ran- 
dom Hamiltonian of the fc-body type. We are interested 
in the way in which properties of the resulting ensemble 
of random matrices affect both, the shape of the average 
spectrum and the spectral fluctuation properties. Our 
approach differs from previous ones in that we are guided 
by an analogy to the metal-insulator transition (MIT). 
In disordered metals, the degree of disorder determines 
whether a system is an insulator (with Poissonian level 
statistics) or a metal (with Wigner-Dyson level statis- 
tics). As disorder increases, a transition from the metallic 
to the insulating regime takes place Q . At the transition 
point, the spectral fluctuation properties are governed 
by a "critical statistic" which usually is characterized by 
three measures 0: (i) The distribution of spacings s of 
nearest eigenvalues is linear in s for small s and falls off 
exponentially for large s; (ii) the variance ^ 2 \L) of the 
number of eigenvalues in an interval of length L is ap- 
proximately logarithmic for small L and linear for large 
L\ (iii) the eigenfunctions display fractional statistics, 
and the distribution of the inverse participation ratios 
(sum of the fourth power of the expansion coefficients 
of the eigenfunctions in an arbitrary basis, see Eq. ^ 
below) scales with N like an inverse fractional power of 
N. Properties (i) and (ii) obviously interpolate between 
Wigner-Dyson statistics for short spectral distances and 
Poisson statistics for larger spectral distances. 

Within the framework of random-matrix theory, an 
ensemble has been constructed that simulates the 
critical statistic at the MIT. That ensemble is a special 
case of a power-law random band matrix (PLRBM). In 
such random matrix ensembles, the variances of the non- 
diagonal matrix elements fall off with some power 
2a of the distance \fj, — u\ from the main diagonal. For 
a > 1 the spectral statistics of the PLRBM is Poissonian, 
for a < I it is of Wigner-Dyson type, and for a = 1 it 
is critical We note that in contrast to Hamiltonian 
matrices for interacting Fermi systems the PLRBM is not 
sparse. 
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Motivated by the analogy to the MIT and by the 
PLRBM in Ref. [loj . we construct a random-matrix en- 
semble (the "scaffolding ensemble" [ScE]) with the fol- 
lowing properties, (i) ScE is more sparse than the Hamil- 
tonian matrix of a fermionic many-body system with k— 
body interactions, (ii) The spectral fluctuation proper- 
ties of ScE are those of the critical ensemble. This then 
suggests (and we show) that for all k > 2, the Hamilto- 
nian matrix of the fermionic problem lies on the metallic 
side of the MIT and is, therefore, chaotic. The case k = 1 
is special and supplies additional arguments that support 
our reasoning. 

Not surprisingly we cannot offer strict analytical proofs 
for some of these statements. Our arguments are based 
on a combination of analytical arguments, numerical ev- 
idence, and the application of a criterion due to Levi- 
tov that is introduced below. 



II. EMBEDDED ENSEMBLE (EGOE(fc)) 

Specifically we investigate a paradigmatic random- 
matrix model that simulates a fermionic many-body sys- 
tem: The Embedded Gaussian Orthogonal Ensemble 
(EGOE) [13]. We consider to spinless Fermions in I > to 
degenerate single-particle states labeled j = 1, . . . ,1 with 
associated creation and annihilation operators a] and , 
respectively. The states carry no further quantum num- 
bers. With 1 < k < to the fc-body Hamiltonian is 



k\ 2 ^ 

3\-3kJ\ J. 
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(1) 



The matrix v is real-symmetric. The matrix elements are 
antisymmetric under the exchange of any pair of primed 
or unprimed indices and are uncorrelated random vari- 
ables with a Gaussian probability distribution with zero 
mean value and a common second moment which with- 
out loss of generality is taken to be unity. The matrix 
representation of in the space of m-body Slater de- 
terminants labeled /i or v with m-body matrix elements 
(u\H^\fj) = H ufl defines an ensemble of real random 
matrices of dimension N = (^) referred to as EGOE(fc). 
We study the spectral properties of EGOE(fc) in the limit 
N ^> 1. For fixed k that limit is reached by letting 
l,m — > oo. Without distinction we consider both the 
dilute limit (m/l — > 0) and the dense limit (m/l — > con- 
stant ^= 0). In Dyson's classification, H^ k ' has orthogonal 
symmetry. Our arguments apply likewise to the cases of 
unitary and symplectic symmetry. As mentioned above, 
for fc < to EGOE(fc) has withstood all attempts at a 
direct analytical treatment [H, [H| . 



III. SCAFFOLDING ENSEMBLE (SCE) 

We first describe the construction of the scaffolding en- 
semble. A criterion due to Levitov then suggests and nu- 



merical simulations confirm that ScE possesses the criti- 
cal behavior characteristic of the MIT. 

We define ScE with the help of a real and symmetric 
scaffolding matrix that has dimension N — 2™ with 
n positive integer. We recall the definition of the auxil- 
iary diagonal for matrices of dimension N: The auxiliary 
diagonal has matrix elements with indices (fx, N + 1 — (J,) 
and \x = 1, . . . , N. We construct the matrix A^ n > by in- 
duction: For n — 1, has dimension 2, zero diagonal 
elements and unit entries in the auxiliary diagonal. Given 
^(n-i)^ -j-jjg ^ WQ diagonal blocks (dimension 2(™ -1 )) of the 
matrix are each occupied by A^ n x \ The elements 
in the two off-diagonal blocks are all zero except for the 
auxiliary diagonal for which all elements have the value 
unity. We display the scaffolding matrix Aw for n = 3 
as an example, 



fO 1 1 1\ 

1 1 1 

10 10 10 

10 10 10 

1 1 1 

10 10 10 

1 1 1 

\l 1 1 0/ 



A (3) 



(2) 



By construction, the matrices A w have the impor- 



tant property ^ w =i ^ 



Y^iv An) 

/i = 1, . . . ,N. Thus, in every row and column of A^ n \ 
the number of non-zero non-diagonal elements is n = 
In N/ In 2. Hence, with increasing distance |/i — v\ from 
the main diagonal, the average density of non-diagonal 
elements of AW falls off like 1/]//- v\. The matrices AW 
share that importa nt p roperty with the power-law ran- 
dom band matrices [lj| that simulate the MIT. Moreover, 
the matrices A^ bear some similarity to the "ultramet- 
ric" matrices studied recently [Til ]. 

With the help of A^ n \ ScE is defined for every n as 
an ensemble ff w of random matrices. The non-zero 
elements of H w reside on the unit elements of the matrix 
A( n ) (hence the name "scaffolding matrix" for A^) and 
on the main diagonal. Except for the symmetry condition 
H$ = H^} the matrix elements are uncorrelated real 
random variables with a Gaussian distribution and zero 
mean values. With a positive and 



n for all 



aS u 



A {n) 



(3) 



the variances are given by 



(1 + 5^)(H$HjZ>) = (S^S vp + 5 MP 5 va )B$ . (4) 

In particular, the variance of every diagonal element is 
equal to a. 

We address the spectral properties of ScE. With E 
the energy and G(E) = l/(E+ - fl») the retarded 
Green's function, the average level density is p(E) = 
— (l/ir)^s(G(E)). Here and in what follows, angular 
brackets denote the ensemble average. To calculate 
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(G(E)), we expand G(E) in powers of and use Wick 
contraction in each term of the sum. Following Ref. [l2[ 
we denote Wick-contracted pairs of matrix elements by 
the same letter and distinguish nested and cross-linked 
contributions. Among the sixth-order contributions, for 
instance, ABCCBA and ABBACC are nested while 
ABC ABC and ABABCC are cross-linked. For n > a, 
only nested contributions contribute to (G(E)). That 
rule (for the case of the GOE demonstrated in Ref. [Hj]) 
can be inferred by comparing the values of low-order 
terms like ABBA, ABBACC, and ABC ABC. Resum- 
mation of the nested contributions gives the Pastur equa- 
tion (G(E)) = (1/E) + (1/E){H^{G(E))H^){G(E)). 
We use Eq. (0} and find 

i E E4p (G(E)) pp ](G(E))^ = V ■ (5) 
p 

To solve Eq. (|5|) we observe that (G(E)} is expected to 
be an analytic function in E with a finite number of 
branch points but without singularity at E — oo. There- 
fore, we expand (G(E)) for \E\ > 1 in a Laurent se- 
ries, (G{E)) pu = Y^Lo Ej ~ P 9$- Inserting that into 
Eq. ([5]) and comparing powers of E we find that non- 
vanishing solutions exist only for j — ±1. For both solu- 

(r>) 

tions we find that the coefficients g^J are proportional 
to the unit matrix for all p. That conclusion hinges 
in an essential way on the fact that Afa) = n for 
all /i so that J2 P ^p = a + n for all fi, see Eq. 
Thus {G{E)) 11U = 5 pv g(E). To determine g(E) we use 
Eq. © and find with A sc = (n + a) 1 ' 2 that \ sc g(E) = 
(E/(2\ sc )) ± - (E/(2X SC )) 2 . The two solutions with 
j = ±1 correspond to the two signs in front of the 
square root. We conclude that the average spectrum has 
the shape of Wigner's semicircle, half the GOE radius 
A GOE cx VN being replaced by A sc = yjn + a w ^fn. 

We have not been able to establish the spectral fluctu- 
ation properties of ScE (and of EGOE(fc)) analytically. 
For the EGOE(fc) we instead use a criterion established 
in Ref. [ll|. To test the applicability of that criterion 
to sparse random matrices, we apply it to ScE. Lev- 
itov investigated a class of random matrices Hfa) for 
which the variances {\H^}\ 2 ) of the diagonal elements 
are much larger than the variances (\Hj^}\ 2 ) (p ^= v) 
of the non-diagonal elements. All of the latter differ 
from zero. Levitov considered the sum S = $^,(1 — 

S^ilH^l 2 ) 1 / 2 / (iH^Jl 2 ) 1 / 2 in the limit of large matrix 
dimension N. Using renormalization-group arguments, 
he distinguished three cases: (i) 5* has a finite limit. 
Then, the spectral statistics of the ensemble is Poisso- 
nian. (ii) S diverges more strongly than In AT. Then, the 
spectral statistics is of Wigner-Dyson type, (iii) S di- 
verges like In A. Then the spectral statistics is that of 
the critical ensemble at the metal-insulator transition. 

ScE and EGOE(fc) are sparse random-matrix ensem- 
bles, and the renormalization-group argument used to es- 
tablish Levitov's criterion is not readily applicable. Thus 



it is not clear whether the criterion applies. If it does, 
ScE must for a 1 ! 2 ^> 1 possess the spectral statistics of 
the critical ensemble since ^ A$ = n = In N/ In 2. In a 
finite system the interesting regime is thus n 3> \fa 3> 1 . 
Fortunately, this regime is just within reach of present- 
day numerical simulations. 

We construct a number of realizations of the ScE 
and compute spectra via matrix diagonalization. We 
are limited to n < 13, and we consider the ScE for 
n = 10, . . . , 13. Our random-matrix ensembles consist 
of about 50 realizations for n = 13 and of up to 30,000 
realizations for n = 10. The wave-function statistics re- 
quires the largest ensembles. For both spectral statistics 
and wave-function statistics, we employ only the levels 
that are in the central 20 % of the spectral densities. In 
this window, the average level density is maximal and 
nearly constant, and finite-size effects can be neglected. 
For the long-range T,^ statistic, we verified that the re- 
sults are insensitive to the degree of the polynomial fit of 
the average level density. 

We first vary a at fixed n = 10 and study the nearest- 
neighbor spacing distribution. Figure [T] shows the tran- 
sition from a delocalized Wigner-Dyson regime to a lo- 
calized Poisson regime as a is increased from a = 3 to 
a = 100. (The ensembles consist of 300 realizations for 
a = 3 and of 3000 realizations for the other values of 
a.) Figure [T] suggests that we are in the critical regime 
somewhere around a ~ 10. 
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FIG. 1: (Color online) Nearest-neighbor spacing distribu- 
tion P(s) of the scaffolding ensemble for n = 10 and a = 
3, 10,30, 100 (histograms from right to left) versus s, the ac- 
tual level spacing in units of the mean level spacing, com- 
pared to the Poisson distribution and Wigner's surmise (i.e., 
the GOE). 

More detailed and extensive numerical calculations 
show that we approach the critical regime (the MIT) for 
a w 17 and n > 13. (Within the available computational 
resources, these parameter values optimally satisfy the 
conditions n 3> \/a 3> 1.) The nearest-neighbor spac- 
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ing distribution shown in Fig. [2] and the long-range E^ 2 ) 
statistic shown in Fig. [3] support this claim. For these 
parameter values, the levels of the scaffolding ensemble 
are clearly strongly correlated at short distances and only 
weakly correlated at longer distances. 



n=13, oc=17 

GOE 

— Poisson 




FIG. 2: (Color online) Nearest-neighbor spacing distribution 
P(s) as in Fig. [T]of the scaffolding ensemble for n = 13 and 
a = 17 compared to the Poisson distribution and Wigner's 
surmise for the GOE. 



Figure [4] shows the distribution of the IPR for the eigen- 
states of the scaffolding ensemble. At fixed a = 17 and 
with increasing n, the distribution becomes scale invari- 
ant (the form of the distribution becomes independent 
of the dimension N = 2™ of the ensemble). Following 
Ref. [HI, we determine the fractal dimension D2 from 
the shift of the IPR distribution that results from a dou- 
bling of the dimension as D2 ~ 0.85. For the largest 
n shown, the IPR distribution exhibits a power-law tail 
with exponent x 2 ~ 1.5. 
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FIG. 3: (Color online) Long-range E^ 2 - 1 statistic of the scaf- 
folding ensemble versus the length L of the energy interval 
(in units of the mean level spacing) for n = 13 and a = 17 
compared to that of a Poissonian spectrum and the GOE. 



FIG. 4: (Color online) Distribution of inverse participation 
ratios (IPR) for a — 17 and n = 10, 11, 12, 13. As n is in- 
creased, the distribution approaches a scale-invariant form. 
The figures behind the symbol "ens" give the number of real- 
izations. 

Let us summarize the main results of this Section. We 
have shown analytically that the average spectrum of 
ScE has the shape of a semicircle. For the three fluc- 
tuation measures that characterize the MIT, our numeri- 
cal results indicate that ScE possesses critical statistic for 
a 1 and n — > 00. This conclusion is on somewhat safer 
grounds than are statements for EGOE(fc) simply be- 
cause ScE is much more sparse: For n = 13, the fraction 
of non-zero matrix elements in every row and column is 
less than 2 x 10~ 3 while for EGOE(fc) with k = 2 and 
matrices of similar dimension, it is about two orders of 
magnitude bigger. For the ScE and n ^> a 1 / 2 3> 1, Lev- 
itov's criterion correctly indicates critical statistic. We 
conclude that Levitov's criterion applies to sparse ran- 
dom matrices. We now use that fact to establish the 
spectral fluctuation properties of EGOE(fc) for k > 1. 



The hallmark of the MIT, however, is a scale-invariant 
distribution of the inverse participation ratio (IPR) [H[ . 
The IPR of a normalized state with expansion coefficients 
-01/ is defined as 

IPR = £hM 4 . (6) 



IV. SPECTRAL PROPERTIES OF EGOE(fc) 

We compare the EGOE(/c) defined in Section ILTl with 
ScE as defined in Section ITLTl We observe that EGOE(fc) 
differs from ScE in three important ways, (i) Count- 
ing shows that for all k = 1, . . . , m the number of non- 
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vanishing non-diagonal matrix elements i/ M „, equal in 
every row and every column, is given by Ylp=i (™) C^,™) 
and, thus, for l,m 3> 1 much larger than In TV w mini. 
Hence, for all k EGOE(fc) is much less sparse than ScE. 
Therefore, we expect EGOE(fc) to lie on the delocalized 
side of the MIT and to possess Wigner-Dyson spectral 
statistics, (ii) The number of fc-body matrix elements 
v 3 -} contributing to a hxed m-body matrix element 
Hp V is, in general, bigger than one. For k = 2, for in- 
stance, the m-body matrix element of two Slater deter- 
minants differing in the occupation numbers of orbitals 
1 and 2 but both with occupied orbitals 3, 4, . . . , (m + 1) 
equals Y?j=3 v lj- I n general, the matrix element con- 
necting two Slater determinants that differ in the occu- 
pation numbers of p single-particle states is the sum of 
\k-p) ^ - t>ody matrix elements. We conclude that the 
variances of the diagonal elements are all equal and given 
by 2 ('?) ; those of the non-vanishing non-diagonal ele- 
ments connecting two Slater determinants that differ in 
the occupation numbers of p single-particle states are 
\k~p)- ft follows that for N — > oo the variances of the 
diagonal elements are much bigger than those of all non- 
diagonal elements. That property is assumed by Levi- 
tov's criterion, (hi) m-body matrix elements occurring 
in different rows and columns may be correlated. For 
k = 2, for instance, the matrix element w|| contributes 
to all m-body matrix elements of pairs of Slater deter- 
minants for which the occupation numbers of orbitals (1 
and 2) and (3 and 4) differ while all other occupation 
numbers agree. Correlations occur only among m-body 
matrix elements in different rows and columns because 
the same A; -body matrix element cannot connect a given 
Slater determinant with two different ones. 

Neglect of Correlations. We hrst consider EGOE(fc) 
under neglect of all correlations, i.e., without property 
(hi), so that fc-body matrix elements appearing in dif- 
ferent locations of the m-body EGOE matrix are as- 
sumed to be uncorrelated Gaussian random variables. 
This assumption greatly increases the number of inde- 
pendent random variables and destroys the connection 
between the resultant random-matrix ensemble and the 
random fc-body Hamiltonian of Eq. (TI| . Under this as- 
sumption we conclude immediately that for N — > oo, 
the spectral properties of EGOE(fc) are for all k the 
same as for GOE. The proof proceeds as for ScE. For 
the average level density we use the fact that — 

VXI^VI 2 ) = Ej=i (T-7)e;)C7") Ureases with N 
much more strongly than the variances 2 (™) of the diag- 
onal elements. That property guarantees that the Pas- 
tur equation holds. The analogue of the ScE relation 
B^p — a + n for all /i also holds: The scaffold- 
ing matrix of EGOE(fc) has the same number of non- 
zero entries in every row and in every column. It fol- 
lows as in Section IIIII that the average spectrum has 
the shape of a semicircle. The radius 2A EGOE ( fc - 1 is 
given by (A™OE (k)? = 2( ™ } + ^ (;) ( ,- m) 



For the spectral fluctuations we use Levitov's crite- 
rion. We avoid infinitely large variances for N — > 
oo by rescaling the energy and all matrix elements of 
EGOE(fc) by the factor (2(™))- 1 / 2 . Then the diago- 
nal elements all have unit variance, and the variances 
(T-p) /POfe)] °^ ^ ne non_ diagonal elements all become 
very small as N — > oo. For all k > 1 the critical sum 

S = E^r(r p P ) 1/2 DC; m )/(2(T)) 1/2 diverges more 
strongly with N — > oo than lniV w mlnZ, and Levitov's 
criterion implies that the spectral statistics is of Wigner- 
Dyson type. 

Influence of Correlations. Correlations among m— 
body matrix elements occurring in different rows and 
columns influence the average level density p(E) and the 
spectral correlations in different ways. For p(E), cor- 
relations cause deviations from the semicircular shape. 
Indeed, the Pastur equation is derived under the assump- 
tion that cross-linked contributions are negligible. That 
assumption fails in the presence of correlations, i.e., when 
{H^Hpa) ^ for {/i, v} 7^ {er, p}. Such correlations, 
nonexistent for k = m, become stronger as k decreases, 
attaining a maximum at k = 1. For p{E) correlations 
cause cross-linked contributions to be as important as 
nested ones. Mon and French [l2j], calculating even mo- 
ments of EGUE(fc) in a basis of Slater determinants and 
using the representation of in Eq. (TTJ), have shown 
that in the dilute limit and for k <C m such contribu- 
tions drive p(E) towards a Gaussian. Thus, the shape of 
the average spectrum of EGOE(fc) is expected to change 
from Gaussian form for k = 1 (where correlations are 
strongest) to semicircular shape for k = m (where corre- 
lations are absent). 

We do not expect correlations between m-body matrix 
elements located in different rows and different columns 
of Hp V to influence the spectral fluctuations of EGOE(fc). 
Wigner-Dyson statistics is a robust property of spectra 
caused by level repulsion. Such repulsion is caused by in- 
dividual matrix elements connecting pairs of close-lying 
levels and is independent of the presence of other corre- 
lated matrix elements. Strong support for this qualitative 
argument comes from the study of EGOE(l). 



V. A SPECIAL CASE: EGOE(l) 

The EGOE(I) is special: The real-symmetric matrix 
Vj, can be diagonalized, and the eigenvalues follow Pois- 
son statistics. The m-body matrix Hp V is then diagonal, 
too, each diagonal element being given by a sum of m 
such eigenvalues. For m 3> I such sums are uncorre- 
lated, and the spectrum is Poissonian. That symmetry 
of EGOE(l) is not obvious in the m-body matrix rep- 
resentation. In excluding such a hidden symmetry for 
k > 2 we appeal to the results of numerical diagonaliza- 
tions. Although done for matrices of small dimensions, 
such calculations should have revealed the existence of a 
symmetry. 
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Because of that special feature, the case k = 1 can 
be used to support some of our arguments and conclu- 
sions very nicely. We compare three ensembles, (i) We 
consider the EGOE(l). As is well known 0, Qj], the 
EGOE(l) has an average spectrum that is (nearly) Gaus- 
sian, and the eigenvalues have Poisson statistics, (ii) We 
consider an ensemble that has the same "scaffolding ma- 
trix" as EGOE(l) but for which all independent m-body 
matrix elements are uncorrelated Gaussian-distributed 
real random variables. In this ensemble, the connection 
with the EGOE(l) is severed, integrability is lost, and 
all correlations present in EGOE(l) are destroyed. Be- 
cause of the lack of correlations between the elements of 
the random matrix, we expect the density p(E) to have 
semicircular shape. Fig. [5] shows that this is indeed the 
case. Because of the loss of integrability we also expect 
Wigner-Dyson statistics for the eigenvalues. This expec- 
tation is confirmed in Fig. [5] (iii) A third ensemble is 
generated by randomly redistributing the single-particle 
matrix elements wj, of the Hamiltonian (fTJ) over the non- 
zero elements of the "scaffolding matrix" of EGOE(l). 
This random exchange of one-body matrix elements also 
destroys the connection of the resulting ensemble with 
the Hamiltonian (fTJ) and, thereby, integrability. How- 
ever, it retains the existence of correlations between ma- 
trix elements. For ensemble (iii) we, therefore, expect a 
(nearly) Gaussian form for the average level density p(E) 
but Wigner-Dyson statistics for the eigenvalues. Our nu- 
merical calculations for I = 12 and m = 6 confirm these 
expectations: Figure[7]shows that p(E) is close to a Gaus- 
sian, and Fig. |S] confirms the Wigner-Dyson statistics for 
the spacing distribution. 

We have devoted particular attention to EGOE(l) not 
only because it is special but also because it is much 
closer to the sparse limit than EGOE(fc) with k > 1. For 
practicable matrix dimensions it is difficult to draw valid 
conclusions about the spectral fluctuation properties of 
EGOE(fc) with k > 2 because the matrices are far from 
the sparse limit. Fortunately, correlations are strongest 
for k = 1. Our conclusions hold, therefore, a fortiori for 
k > 2. In particular, the results for EGOE(l) strongly 
support the conclusions drawn at the end of Section IIVI 



VI. CONCLUSIONS 

We have constructed a random-matrix ensemble (ScE) 
that is more sparse than EGOE(fc) for all k. We have 
shown that ScE mimics the metal-insulator transition 
and possesses critical spectral statistics. Using ScE as a 
test case, we have verified that Levitov's criterion applies 
to sparse random matrices. 

Comparison with ScE suggests that for all k > 1, 
EGOE(fc) is on the delocalized side of the metal-insulator 
transition and possesses Wigner-Dyson spectral statis- 
tics. We have presented a number of arguments which 
strongly support that expectation: Levitov's criterion in- 
dicates chaos, and properties of the (modified) EGOE(l) 




FIG. 5: (Color online) Level density p(E) (compared to a 
best-fit semicircular density) for a random-matrix ensemble 
with a scaffolding matrix of the EGOE(l) and random, un- 
correlated matrix elements (case (ii) discussed in the text). 
The ensemble consists of 1000 realizations for I = 12 orbitals 
and m = 6 fermions. 
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FIG. 6: (Color online) Nearest -neighbor spacing distribution 
P(s) compared to Wigner's surmise for a random-matrix en- 
semble with a scaffolding matrix of the EGOE(l) and random, 
uncorrelated matrix elements (case (ii) discussed in the text). 
The ensemble consists of 1000 realizations for I = 12 orbitals 
and m = 6 fermions. 

illustrate the different roles played by integrability of the 
underlying Hamiltonian on the one hand, and the ex- 
istence of correlations between matrix elements on the 
other. We conclude that spectra in finite many-body 
systems governed by few-body interactions generically 
display Wigner-Dyson level statistics. By implication we 
conclude that in the limit of infinite matrix dimension, 
the distribution of the eigenfunctions is Gaussian. These 
conclusions hold for all three symmetry classes (orthog- 
onal, unitary, symplectic). 
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EGOE(fc) is based on the assumption that the single- 
particle states are degenerate. If that degeneracy is lifted, 
chaos may be reduced. We have not specifically ad- 
dressed that case Q. 



FIG. 7: (Color online) Average level density p(E) (compared 
to a best-fit Gaussian density) for a random-matrix ensemble 
with a scaffolding matrix of the EGOE(f ) but with random- 
ized correlated matrix elements from the EGOE(f ) (case (iii) 
discussed in the text). The ensemble consists of 1000 realiza- 
tions for I — 12 orbitals and m — 6 fermions. 




FIG. 8: (Color online) Nearest-neighbor spacing distribution 
P(s) compared to Wigner's surmise for a random-matrix en- 
semble with a scaffolding matrix of the EGOE(l) but with 
randomized correlated matrix elements from the EGOE(l) 
(case (iii) discussed in the text). The ensemble consists of 
1000 realizations for I = 12 orbitals and m = 6 fermions. 



We believe that the results of this paper, although not 
entirely based on strict analytical arguments, convey new 
insight into the mechanisms that determine the spectral 
shape and the spectral fluctuation properties of fcrmionic 
many-body systems. 
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